z = [5;-25;-15;10;-20;-10;25;-5;5];
A = [-1,1,1,0,0,0,0,0,0;-1,3,-1,0,0,0,0,0,0;0,0,0,-3,1,1,0,0,0;0,0,0,-1,1,-1,0,0,0;1,0,0,1,0,0,1,0,0;0,1,0,0,1,0,0,1,0;0,0,1,0,0,1,0,0,1];
b = [0;0;0;0;100;120;80];
Ib = [0;0;0;0;0;0;0;0;0];
[x,fval] = linprog(z,A,b,[],[],Ib)